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Or 



O 

00 



i 



Hiroshi Kuratsuji and Shouhei Kakigi 
Department of Physics, Ritsumeikan University-BKC , Kusatsu City 525-8577, Japan 

(September 18, 1998) 

qq'i A theoretical study is given of a new type of optical vortex in nonlinear anisotropic media. This 

' is called an "optical spin vortex" , which is realized as a special solution of the two-component non- 

| linear Schrodinger equation. The existence of spin vortex is inherent in the spin texture that is 

caused by the anisotropy of dielectric tensor, where a role of spin is played by the Stokes vector. 
By introducing the effective Lagrangian governing the non-linear Schrodinger equation, we derive 
the evolution equation for the motion of spin vortex, which exhibits a clear distinction from the 
• conventional singular vortex inherent in the phase function that is described by the single component 

complex field. 

\ PACS number: 42.25.Ja, 42.65.-k, 67.57.Fg, 78.20.Ek 

The nonlinear optics has been a major subject in physics for long time [Q. The main interest is focussed on 
^ ■ particular solutions of the basic non-linear equation that governs the light field or eletro-magnetic field. Among 
others, remarkable is the optical vortex, the existence of which has been early suggested in |Q. Recently the detailed 
study has been carried out from both of theoretical and experimental point of view |^,|3| . More recently experimental 
verification has been also given for the multi- vortices [Q . The basic idea of the optical vortex follows an analogy with 
the existence of the superfluid vortex that is described by the complex Landau-Ginzburg order parameter for bose 
fluid as is stated in ||. The equation for the light field, which is known as the nonlinear Schrodinger (NLS) equation, 
is very similar to the Pitaevki equation. Thus it is natural to expect the occurrence of the optical counterpart of the 
I | superfluid vortex. 

The purpose of this letter is to put forward a possible new type of vortex and its dynamical behavior in the 
nonlinear and anisotropic media. We call this new type vortex the "optical spin vortex" . Our starting equation is 
the two-component NLS equation. Here the basic ingredient is the state of polarization, which is realized by the 
!— 1 1 two-component field. The concept of polarization is fruitful which forms a basis of modern crystal optics [^|^]- The 
quantity that describes polarization state is given by the Stokes parameters (or vector), which is geometrically realized 
as a point on the Poincare sphere. The crucial point is that the Stokes vectors form the field of pseudo-spin, namely, 
■ we consider the special configuration of pseudo-spin. Having recognized the role of the Stokes vector as a pseudo- 
l/"") [ spin, we have the analogy with the spin vortex. This can be naturally formulated by the effective Lagrangian for 
the pseudo-spin field reduced from two-component NLS equation. Such a Lagrangian has not been used thus far in 
condensed matter physics. Notable here is that the spin vortex is non-singular in contrast to the singular vortex that 
appears in the conventional optical vortex, where the singularity is inherent in the phase of the complex scalar field. 
By using the effective Lagrangian, we can naturally derive the evolution equation for the non-singular vortex, which 
reveals the topological invariant that is characteristics of the non-singular vortex. In this note, we are concerned with 
the basic qualitative idea and the details of numerical analysis will be given elsewhere. 

Two component nonlinear Schrodinger equation. — First we shall derive the two-component nonlinear Schrodinger 
equation for the light wave traveling through anisotropic media. The procedure follows the one developed in a recent 
paper . Suppose that the electro-magnetic wave of wave vector k travels in the direction of z with the dielectric 
tensor e. The nonlinear nature of media implies that e has a field dependence in nonlinear form. The z-axis is 
prescribed to be the principal axis of the dielectric tensor, namely, the axis corresponding to one of the eigenvalues of 
the dielectric tensor. In this geometry, e is taken to be 2 x 2 matrix and the Maxwell equation for the displacement 
field D is given by the second order differential equation: 
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where V = {-§^, -§^) and (x, y) denotes the coordinate in the plane perpendicular to z axis. Now we put 

D(z) = f(x, y, z) enp[iknoz] (2) 

with k = — , and no(= v^o) means the refractive index for the otherwise isotropic medium. The amplitude f(x, y, z) 
is written as f = '(/lj/a) = fi^i + $2^2, which is slowly varying function of z besides the (x,y) dependence, and 
ei and e2 denotes the basis of linear polarization. By subsituting (0) into (111) and considering the short wave length 
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limit, i.e., |^|| << k : we can derive the equation for the amplitude f; namely, keeping only the first derivative ^| 
besides the Laplacian with respect to (x,y), we get 
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where A is the wavelength divided by 2n. This equation is regarded as a two-state Schrodinger equation where A just 
plays a role of the Planck constant and z is a substitute of the time variable. The components (/i,/2) couple each 
other to give rise to the change of polarization which is just the effect of birefringence governed by a 2 x 2 matrix 
"potential" v — e — n\. The v represents the deviation from the isotropic value and if the non-abosortive media is 
considered, it becomes hermitian. From the hermiticity, the most general form of v is written as 



v a + a f3 + ij 
(3 — «7 vq — a 
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Now for the use of later argument, it is convenient to transform the basis of linear polarization f in the circular basis, 
that is, 

e± = -^=(ei ±ie 2 ) (5) 
which defines the unitary transformation (e + ,e_) = T(ei,e2) Here T is the unitary transformation of 2 x 2 matrix: 

Hence we can introduce the wave function and as i/j = Tf = * (ip* , ) . Thus we have 
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dz 

and the transformed "Hamiltonian" becomes 

A 2 

H = ThT- 1 = V 2 + V (8) 

n 

where the potential term is written in terms of the Pauli spin; V = v x 1 + J^j = i UjUj. We here introduce the 
"quantum" Lagrangian leading to the Schrodinger type equation, which is given by 

I = J V f MA-^ - Hj *pd 2 xdz. (9) 

Indeed, the Dirac variation equation 81 = recovers the Scrodinger equation. For later convenience, we write 
L C = J ipUX-^ixfx and H = J^H^cfx, which are called the canonical term and the Hamiltonian term respectively. 

Reduction to the pseudo-spin field. — Now, having given the two component field function ip, we can introduce the 
Stokes parameter by using the Pauli spin J^], namely, this is written as Si = ip^criip, Sq — ip^lip with i — x,y, z . We 
see that the relation Sq = + Sy + S 2 holds. So gives the field strength; Sq = |D| 2 . Using the spinor representation, 

9 9 

V'i = So cos-, i>2 = Sosin-exp^], (10) 

we have the polar form for the Stokes vector S = (S x , S y , S z ) ; 

S x = So sin^cos^, S y = So sin 6 sin <j>, S z — Socos9 (11) 

which forms a pseudo spin, so to speak, and is pictorially given by the point on the Poincare sphere. In terms of the 
angle variable, the canonical and Hamiltonian term in the Lagrangian are given as 

Lc^J^il-cos^d^ 

H = H T + V. (12) 
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Here the potential term in the Hamiltonian is given as 
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V= I '(vo + J^vA^x. (13) 
J t=i 

We note that vq and Uj's are nonlinear functions of the field strength 5*0 in general and this fact suggests that we 
have a stable special solution for the field of pseudo spin. The kinetic energy term is written as a sum of three terms; 
Ht = ^- / Vijj^ipcPx = Hi + i?2 + H3 where the respective terms become 

Hi = / ^(VS ) 2 d 2 x, 
J n 

H 2 = /^{(W) 2 +sin 2 tf(V0) 2 }d 2 x, 
J 4n 

H 3 = f ^-{(1 - cos 9)V<f)} 2 d 2 x (14) 

where Hi gives the kinetic energy associating with the space modulation of the field strength, and the second term 
represents the intrinsic energy for pseudo-spin which exactly coincides with the continuous Heisenberg spin chain [^| . 
The last term is nothing but the fluid kinetic energy inherent spin structure, where we put 

v = (1 - cos6»)V0 (15) 

which means the velocity field as is seen below. In this way, we have a new type of fluid, so to speak, "anisotropic 
optical fluid" , which is naturally described by the Lagrangian written in terms of the pseudo-spin field. 

Equation of motion for optical spin vortex. — We now examine the vortex solution by starting with the pseudo- 
spin Lagrangian given in the above. Taking variation leads to the coupled equations for the field variables (Sq,9, <p). 
In what follows, we confine our argument to the case that Sq becomes constant. Physically, this corresponds to the 
constant background field with a "dark" core which is controlled by the profile of the remaining variables (0,0). A 
static solution for one vortex is obtained with the phase function 4> — tan -1 ^, and the profile function 6 is given as 
a function of the radial variable r. The explicit form of the profile function 0(r) may be derived from the extremum 
condition for the Hamiltonian H, with the specific boundary condition at r = 00 and r = 0. Namely, we choose 
such that at the origin r = 0, we can take 9(0) = in general, that is, the left-handed circular polarization at the 
origin [see Fig.|l|(a) and (b)], whereas at r — 00, there exist several possibilities; we here consider two typical cases: 
A) 9(oo) — 7r and B) 9(oo) — ir/2. We here introduce the vector m(x) = S/Sq, Then, we have 7773(0) = 1 for both 
cases, A), B), namely, the spin field directs upward. On the other hand, we have at r = 00, 7713(00) = —1 for case A) 
and 7713(00) = for case B), that is, the spin field directs downward for A), which represents the right-handed circular 
polarization [Fig.[l](a)] and outward for B), which represents the linear polarization [Fig.|l](b)] . Here we simply assume 
the existence of the solution #(r)satisfying the above boundary condition. 

In order to treat the motion for a single vortex, we introduce the coordinate of the center of vortex, R(z) = 
(X(z),Y(z)), by which the vortex solution is parameterized such that 9(x — R(z)) and </>(x — H(z)). By using the 
parameterization prescribed in the above, the canonical term Lc, the first term in (113) , is written as 



L c = I -^(l-cos6»)V</)-Rd 2 a; 

S ° X I « Rd 2 x (16) 

where we have used the relation: |j = J^;R, = —V(f> with R = 4^. Remarkable point here is that the velocity 
field (|l5| ) does not bear singularity at the origin compared with the behavior of the conventional vortex for which the 
velocity field is simply defined by the gradient of the phase v = V<j) ||. Keeping this in mind, we shall derive the 
equation of motion for a single vortex. In order to achieve this we consider the Euler-Lagrange equation for R, which 
gives the "balance of forces" 

_ d 3L C dL c _ dH 
c dz dR <9R OR' [ ' 

By using eq.(|l5|), we get 
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i-(kxR) = - M (18) 



where the k in ( |l8| ) is the unit vector perpendicular to the xy-plane. Here a is defined as 

< 19 > 

In deriving ( jig) , we have used the relation = — . The integrand of a is nothing but the vorticity which we put 
u>. Using eq.(|15|), we can write u> in terms of the angular variable: 

u = V x v = sin 9(V8 x V</>) (20) 

which is also rewritten in terms of the spin field m 

VXV - m '(^ X ^j- (21) 

Equation (|l]) (or (|2p|)) is an optical counterpart of a topological invariant of hydrodynamical origin jll]] or otherwise 
the Mermin-Ho relation for supcrfluid Hc3- A |L2] . We get an integral of the differential form; 



sin Odd A d<f>. (22) 
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Here taking into account of the boundary conditions for the profile function 8(r) mentioned above, we have a topological 
interpretation for a. In the case A), the vortex configuration gives the continuous mapping from the compactified 
space Si to the spin configuration S2, so the a in has a meaning of the mapping degree of S2 — > Si. Hence 
we get the quantization of a: a = to (m=integer) . For the case B), on the other hand, the mapping becomes 
S2 — > Si/2 (hemisphere), so intuitively we have the quantization a = to/2 (to: integer). Thus it should be noted 
that this distinguishes the topological invariant for the conventional optical vortex [Bj, which gives the invariant of 
§ c S7<j>ds = 2-7T. It should be noted that the left hand side of the equation of motion ( |l8| ) is considered to be an 
optical analog of the Magnus force in superfluids or superconductors. As is mentioned above, the "force" incorporates 
the topological structure inherent in the coreless vortex, which is quite different from the conventional optical vortex 
which is singular at the origin. The point here is that there are two distinct types vortices according to the different 
boundary conditions at infinity. 

The difference between two types of vortex is significant from the point of view of condensed matter physics, namely, 
we expect that the non-singular spin vortex is converted to the usual phase vortex. This may be considered to be a 
kind of phase transition as a consequence of the phase transition from anisotropic and isotropic that occurs in the 
medium. Actually, as a result of transition the birefringence disappears and the Hamiltonian (or dielectric tensor) 
becomes scalar matrix, which means that the two component equation is converted to one component; ip = Sq exp[i<j>], 
where the amplitude vanishes at the vortex center. In the process of phase transition the core structure of the optical 
vortex may survive and it may be possible to observe the structure change between two cores; the soft and singular 
[SeeFig.§. 

The case of multi-vortices. — We now consider the correlated motion of two or more vortices. The canonical term 
in Lagrangian is simply obtained by replacement as v • R — > . • Ri where Vj is the velocity field coming from the 
i-th vortex. The effective Hamiltonian for the assembly of vortices is given by the fluid kinetic energy, Jl3| namely, 

Heff= /5>(y')Vlog|y'-(*-Ri)l 

3 ij 

x ^(y")Vlog \y" - (x - R J )\dP V , dPy"dPx. (23) 



We can intuitively guess that it depends on the distance between center coordinate Ry = R s ; — Hj. The effective 
Lagrangian is written as 

T S*\ ^( dXi dYA 

i v / 

Here a special case, we consider the equation of motion for a couple of vortices, which becomes 
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Sq A _ dRi\ dH e ff 



2"l kx ^rJ^^7 (25) 

with i — 1,2. By introducing the relative coordinate r = Hi — R 2 ,we have a reduced form 

^(kx^U-^. (26) 
2 \ dz) dr K ' 

From this we obtain the constant of motion r 2 — C , which is reminiscent of the well known conservation law of 
the vortex dynamics. Furthermore, we have the relation between the polar angle and the relative vector r, namely, 
d> = f(r) = constant, which shows that the two vortices rotate around each other by constant angular velocity (c.f. 
) , [Fig.^l . Finally, we give a remark on a possible conversion of the effective Hamiltonian in accordance with the 
change between two types vortices (non-singular and singular types). Namely, for the singular vortex it is given by 
the logarithmic form H e ff = — ^p.^ log |Rj — Rj|. On the other hand, in the limiting case that |y'| << |x — Rj|, 
the effective Hamiltonian ( p3| ) is shown to be reduced to the form 

Vlog|x|Vlog|x- R l3 \d 2 x ~ cr 2 log|Ry-| (27) 

which implies that the Hamiltnian has the same form as the singular case except the coupling constant in such a 
limiting case. The difference of the coupling constant may play a substantial role for detecting the possible structure 
change between two types of optical vortices. 
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FIG. 1. Profile of the non-singular vortex; a) The case 
of 0(oc) = ty and b) The case of 0(oo) = f . 



FIG. 2. Illustration of for evolution of vortex pairs and 
the possible transition between the non-singular (b) and 
the singular (a) vortex. The bullet • and white circle O i n 
the right hand side indicate an instantaneous vortex pair 
at a particular value of z. 
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